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. Abstract 

Acoustic black holes are very interesting non-gravitational objects which can be described 
""q I by the geometrical formalism of General Relativity. These models can be useful to ex- 

perimentally test effects otherwise undetectable, as for example the Hawking radiation. 
^ { | The back-reaction effects on the background quantities induced by the analogue Hawking 

bi), radiation could be the key to indirectly observe it. 

We briefly show how this analogy works and derive the backreaction equations for the lin- 
earized quantum fluctuations in the background of an acoustic black hole. A first order in 
h solution is given in the near horizon region. It indicates that acoustic black holes, unlike 
Schwarzschild ones, get cooler as they radiate phonons. They show remarkable analogies 
with near-extremal Reissner-Nordstrom black holes. 
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1 Introduction 



Everybody agrees on the fundamental importance in Physics of experimental 
confermations of the theoretical predictions. Yet unfortunately some fields 
of Physics suffer an incredible lack of experimental data, which can be a 
serious obstacle for their future developments. Among these is Quantum 
Field Theory in curved space (the study of quantum fields in presence of 
gravitational sources), whose most spectacular prediction - the Hawking ef- 
fect - widely considered a milestome in theoretical Physics, so far has no 
experimental evidences at all which could confirm its existence. Taking into 
account Quantum Mechanics, Hawking [1. showed that black holes are unsta- 
ble objects which thermally radiate particles at a temperature proportional 
to the horizon surface gravity k, i.e. Th = 2 ^k b c * wnere c * s the speed of 
light and Kb the Boltzmann constant. Because of its weakness, it is rather 
impossible to detect Hawking radiation from astrophysical observations: the 
emission temperature for a solar mass black hole would be ~ 10 -8 K, far 
beyond the cosmic microwave background temperature of 2.7 K. Moreover 
things become even worst if the mass of the hole is higher. On the other 
hand, also detection of Hawking radiation from primordial mini black holes 
seems problematic since inflation would have dilute them together with all 
the population of primordial objects [2j. 

A key feature of Hawking radiation is its absolute independence from the 
dynamical details. In fact it is not directly linked to the Einstein equations, 
but only to the geometrical properties of space-time near the horizon. It is 
a pure kinematical effect, so in principle it will occur whenever a space-time 
structure similar to that of a black hole appears. 

So recently a gleam has been opened in indirect detecting Hawking radia- 
tion through the so called analogue models of gravitational systems in Con- 
densed Matter Physics jHJSE]- It seems possible to simulate in laboratory 
black hole-like objects thanks to the astonishing analogy which exists be- 
tween sound propagation in an inhomogeneous fluid flow and light trapped 
in the gravitational field of a black hole. In the following we will focus the 
attention on a hydrodynamical flow in a nozzle, which is the simplest setting 
where the analogy can work. In this kind of systems, it should be possible to 
create in laboratory acoustic analogue of gravitational black holes which shall 

$A hint of its fundamental origin can be given observing this formula: three universal 
constants (h from Quantum Mechanics, c from Special Relativity and G, contained into 
k, from General Relativity), appear. 
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thermally radiate phonons for the same reason that gravitational black holes 
emit Hawking radiation, providing a new context to detect, albeit indirectly, 
it. 



2 The analogy 

Let us consider a fluid which is irrotational (i.e. the velocity v can be ex- 
pressed as the gradient of a potential ■0) and homentropic (i.e. the pressure P 
is a function of the density p only). The equations of motions which describe 
its flow are the continuity and the Bernoulli equations 

p + V-(pv)=0, (1) 
iP + ±i? + p(p) = J (2) 



where p(p) = with u(p) the internal energy. 

Now, expanding the fields ip, p, v around a background solution of Eqs. (JUEJ) 
as ip = ipQ + eipi, p = po + epi, v = v + evi up to 0(e 2 ), with e a dimensionless 
expansion parameter, the first order in e equations of motions read 



Pi + V ■ (piVq + pqVx) = , (3) 
ipi + v ■ v x + Jpi = , (4) 



where we have introduced the speed of sound c as c 2 = p^\ Po - Inserting pi, 
extracted from Eq. (JU), into Eq. Q it yields 

-d t [^(dt^+vo- V^i)]+V-{^ +p„VVi} = 0. 

(5) 

Since sound is the first order fluctuation of the quantities which describe the 
mean flow, Eqs. (JU) and (JSJ) describe the usual propagation of sound waves. 
The crucial point is that Eq. © can be written as the four-dimensional 
d'Alambertian equation 

fy(V=^^i) = 0, (6) 

provided one introduces an "acoustic metric tensor" 

= Pof -(c 2 - vl) -i% \ m 
9flu ~ c \ -vo 1 ' 1 ] 
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Figure 1: Schematic representation of a de Laval nozzle. The fluid is supposed 
to flow from right to left, reaching the sound speed c at the waist. 

So the propagation of sound waves coincides with the propagation of a mass- 
less scalar field in a fictitious curved background given by the acoustic met- 
ric g^ v . It is worth noting that this is not the real Minkowski space-time 
metric lying underneath the fluid flow, which is completely non-relativistic 
(| «o| <S cught)- g^v is the metric the fluctuations "feel" and is due to the non- 
homogeneities of the fluid flow which, as it should be, modify and distorce 
the motion of sound waves 

Inspection of (|ZJ) reveals that the component goo vanishes for \v\ = c. This 
fact suggests that not only a four-dimensional Lorentzian formulation of the 
sound propagation can be given, but also that this propagation, when the 
flow reaches the speed of sound, can be similar to that of light in a black hole 
space-time. 

The simplest device to built this " acoustic black hole" is a de Laval nozzle , 
a converging- diverging nozzle (see FigHJ): let us suppose that the fluid flows 
from right to left. The flow can be adjusted to be subsonic in the right region, 
it accelerates reaching the sound velocity exactly at the waist, then becomes 
supersonic into the left [7] . It is clear that for upstream sound waves propa- 
gating in the supersonic region of a de Laval nozzle is impossible to reach the 
subsonic region on the right of the nozzle. They are dragged by the flow and 
behave as light trapped in the gravitational field of a black hole. Therefore 
the supersonic side of the nozzle represents the sonic black hole and the waist 
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the sonic horizon. In addition it is possible to show that the sound moves 
following the geodesies of the acoustic black hole metric (J7J). Now, pushing 
forward the analogy, after quantization of the field ip\ and using the same 
Hawking arguments, a flux of phonons (the sonic quanta of the field ^i) at 
the temperature Th is expected: 

T» = ^- (8) 



where here c is the speed of sound and k is the surface gravity calculated on 
the acoustic horizon defined as k = — | d<yC ^ ^ [HI IE] This is the sonic 
analogue of Hawking radiation. In some cases this temperature is not so far 
from the actual experimental limits (see |HE])- 



3 Back-reaction of phonons 

Unfortunately this flux of phonons, because of the extremely low associated 
emitted power (~ 10~ 31 W for water, ~ 10 -28 W for Bose-Einstein conden- 
sates), is still rather difficult to display; so it is of great interest to know how 
the phonon flux influences the background flow, i.e. to know its back-reaction 
effects, in order to have a chance to detect indirectly the Hawking radiation 
by the observation of its back-effects. But this is not a trivial extension of the 
analogy, since the back-reaction involves the dynamics of the system, while 
the analogy works only at the kinematical level. This issue can be addressed 
by considering the expansion of the equations of motion up to second order 
(i.e. taking ip = ip + eipi + e 2 ^, p = po + epi + e 2 p 2 , v = v + evi + e 2 v 2 ). 
The first order equations of motion (0 EJ) contain linear fluctuations (they 
describe the propagation of sound waves), while the back-reaction is given 
by the second order equations, where the fluctuations appear quadratically: 

h + V • (p 2 v + pqv 2 ) + V • {pivi) = , (9) 

1 c 2 

^2 + vo-v 2 + p(p2) + -v 2 1 -—p 2 1 = , (10) 

L Zpo 

where we have assumed for simplicity sake c = const. 

Making use of the standard definition of the stress-energy tensor^ for a mass- 
less scalar field T M „ = d^idytpi - \gf J ,ud a ip 1 d ot ijji and Eq. (gj), Eqs. (P fTU|) 

§To be precise is not the stress-energy tensor which gives the energy-momentum 
of the whole flux. It just refers to the energy-momentum content due to the fluctuations 



4 



read (after quantum averaging over the fluctuations): 



h + V • (p 2 v + p v 2 ) + V • 

^2+V -V 2 + /i(p 2 ) - ~ ('- 



Pb 



Wii)) =0, (11) 



c 



(T) = . 



(12) 



where T = g al3 T a p is the trace of the pseudo stress-energy tensor and ( ) must 
be intended as quantum expectation values taken on the quantum state which 
properly describes the formation of a black hole, i.e. the Unruh state [TT] . 
It is worth noting that, although the dynamics for an hydrodynamical fluid 
flow deeply differs from the gravitational one, back-reaction effects are driven 
in both cases by the stress energy tensor, even if its components appear in 
different ways. Continuity equation, as one could have expected, acquires an 
extra term which takes into account the momentum-density of the phonons 
(—^r(T ti + v J B Tij) = p\v\ = yf^gT®). The Bernoulli equation instead, which 
provides the modification to the velocity, gets modified by the trace of the 
stress energy tensor. In the gravitational analogue equation is the flux which 
determines the mass correction. This is not surprising since one can not 
expect the very same back-reaction equations for the two settings, that, even 
if identical from the kinematical point of view, are completely different from 
the dynamical one [T2"l ITS} ITIj. 

4 An approximated solution 

A general solution for Eqs. (fTH [T2"|) can not be given in a general four- 
dimensional setting, since the explicit form of the stress-energy tensor com- 
ponents is unknown. The problem can be partially addressed by introducing 
the approximation of quasi-one dimensional flow (i.e. the transverse - let say 
x, y - components of the velocity are negligible with respect to the component 
along the axes of the nozzle - say z), which allows to perform a dimensional 
reduction. So the two-dimensional back-reaction equations become 



(i.e. the sound waves). For this reason it is usually called "pseudo" stress-energy tensor. 
For a detailed treatment see |H1 IIP). 



Ap 2 + d z [A(p V 2 + p 2 V )] + d. 



h\-- c ((T t ^) + v (T^)) 



0, (13) 



A(ij 2 + v -v 2 + p(p 2 ^--(T^) = 



(14) 
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Here A(z) is the area of the transverse section of the nozzle and T 



(2) 



(lb 



A(d a ipidbipi — ^g^ cd d c ipidd'ipig ab ), with g ab } the (i, z)-section of the four di- 
mensional acoustic metric (0). 

The quantum field ip% now satisfies a two-dimensional equation for a mass- 
less scalar field coupled not only with the two-dimensional metric but also 
to the combination Ap/c, which reminds the four-dimensional origin of the 

model: d a \/— g^g ab dbipi J = 0. This coupling causes backscattering 
of the modes. Now, neglecting this backscattering in accordance with the 
Polyakov approximation ^H] , the components of the stress-energy tensor are 
known for any general two-dimensional metric (see [Tfi]). 
We are interested in the first order in h correction to the mean flow, so 
the (Tab) components must be evaluated on the background flow. More- 
over we have chosen a shape for the transverse area of the nozzle as A(z) = 
Ah{1 + k, 2 z 2 ), which is always valid if one restricts enough his analysis in a 
proper near horizon region. Ajj is the area at the waist, posed in z = 0. For 
this particular choice for A, the horizon surface gravity becomes k = c 2 k. 
The quantum corrections to the velocity and density are: 

^ = A%w-Wc ^ + * Kz)Kt > (15) 
g2p2 = A\ple-^c PHe c 2 {a + 6KZ)t ' (16) 



where b\ = 97/2, c\ = —3047/15, a = —a 2 + cnbi, 5 = —b 2 — a 2 + ckci, 
with a 2 = 11897ck/120, b 2 = -1618177ck/120 and 7 = Kc 2 A^/24vr. These 
solutions are only valid in the near horizon region (kz -C 1) and for small 
times (cut 1). That means that they just give the hints of how the back- 
reaction starts modifying the background flow near the horizon. Inspection 
of Eqs. ([TBI ITBj) reveals a depletion of both velocity and density. 
Then the position of the horizon, given by \v correctec i\ = c, moves: z H = 
— — < 0, shrinking the sonic black hole. 
Also the emission temperature starts changing: 



m h dv 
T 



2,tc k r dz 



he 

; 

2txkb 



563 * A ' 

1 €K cA H t 

720vr 



(17) 

4TYKB f^U7i 

As time goes on, the temperature of the emitted phonons decreases. Extrap- 
olating this result supposing the system passing through stationary quasi- 
equilibrium states, one can try to get the late time evolution of T H . What 
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we have found is a lowering of T# oc t~ 2 : the sonic black hole tempera- 
ture vanishes in an infinite amount of time. This behavior completely differs 
from what one finds in Schwarzschild black holes, where the back-reaction 
increases the temperature and the black hole lifetime is of order Mq 3 (Mq is 
the initial mass of the hole). Acoustic black holes instead show a behavior 
similar to charged Reissner-Nordstrom black holes, whose temperature tends 
to zero as oc t" 3 [THllli] . 

5 Comments 

The acoustic analogy in Condensed Matter Physics of gravitational objects is 
probably one of the most hopeful way to understand and verify the Quantum 
Field Theory (QFT) in black holes space-time (in particular Hawking radia- 
tion). Even if the first article by Unruh [3| is more than twenty years old, a 
lot of work has still to be done. In fact, no experimental proofs of Hawking 
radiation has been found yet, since the effects to point out are still quite 
tiny, even if not completely out of the present (or near future) experimental 
capability. 

An easier way to experimentally observe these kind of effects could be to ap- 
preciate the indirect presence of them via the modifications they induce on 
the underlying system. So the study of back-reaction could be of primarily 
importance. In this work we have briefly presented the very first attempt in 
the analysis of back-reaction in the simplest acoustic setting that provides a 
black hole analogue: a supersonic fluid flow in a de Laval nozzle. It displays 
that if the dynamics which drives the analogue system is correctly described 
by hydrodynamical equations, the back-reaction due to the Hawking radia- 
tion is governed by the stress-energy tensor of the phonons.^ 
A solution, under some quite reasonable approximations, tells us that the 
horizon should move shrinking in size the sonic black hole and that the tem- 
perature emission should drop till to completely vanish in an infinite amount 
of time. 

This is only a preliminary conclusion, that has to be confirmed by numer- 
ical analysis and by further deeper studies, becoming the first step in the 
direction to fully make out the real prospects contained into the analogy. 

^This is not a general result: it holds only if the dynamics is properly described by the 
hydrodynamical equations. For example for Bosc-Einstein condensates it will not be true 

en 
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